Scale Separation Scheme for Simulating Superfluid Turbulence: Kelvin wave Cascade 
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A Kolmogorov-type cascade of Kelvin waves — the distortion waves on vortex lines — plays a key 
part in the relaxation of superfluid turbulence at low temperatures. We propose an efflcient numeric 
scheme for simulating the Kelvin wave cascade on a single vortex line. The idea is likely to be 
generalizable for a full-scale simulation of different regimes of superfluid turbulence. With the new 
scheme, we are able to unambiguously resolve the cascade spectrum exponent, and thus to settle 
the controversy between recent simulations Q and recently developed analytic theory 0- 



PACS numbers: 67.40.Vs, 03.75.Lm, 47.32. Cc 



The superfluid turbulence (for an introduction see 
0, 0) can be viewed as a tangle of quantized vortex 
lines. In the case of an essentially finite temperature, the 
superfluid turbulence is perfectly described by the local 
induction approximation (LI A), as it has been shown in 
pioneering simulations by Schwarz 0. LIA reduces 
the problem of finding the velocity of an element of the 
vortex line to its local differential characteristics: 



s = /3 s' X s" 



friction 



/? = — HR/ao) 
in 
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where s = s(^, t) is the evolving in time radius-vector 
of the vortex line element, the parameter ^ being the 
arc length, the dot and prime denote differentiation with 
respect to time and ^, respectively; k is the quantum 
of velocity circulation; R is the typical curvature radius 
(which is of the same order that the typical interline spac- 
ing), and ao is the vortex core radius. A necessary con- 
dition of applicability of LIA is 



HR/ao) » 1 



(3) 



At zero temperature, a crucial role in the dynamics 
of superfluid turbulence is played by the Kelvin wave 
cascade 0, S B B 03 • LIA iFails to properly describe 
a considerable part of the cascade inertial range (see, 
e.g., and references therein). In this case, one has to 
use the description in terms of full Biot-Savart equation 
(BSE): 



s) X dso/|so — s|' 



(4) 



Here the evolving in time radius- vector of the vortex line 
element, s = s{(,t), is labelled by some parameter C 
(not necessarily the arc length), the vector sq is of the 
same physical meaning as s, understood as an integra- 
tion variable; the integration is over all the vortex lines; 
dso = (dso/dC) d(^- The two conditions of applicability 
of Eq. are: (i) the absence of the normal component 
implying the absence of the friction term, and (ii) 



i? > ao 



(5) 



The need to use the full Biot-Savart law in dealing 
with the Kelvin wave cascade had been recognized in 
Refs. It is worth noting that the Kelvin wave cas- 

cade is not the only problem where non-local vortex line 
interactions play a crucial part. The other examples in- 
clude the conditions of stability of vortex knots revealed 
by Ricca et al. [131: the vortex knots that are unstable 
under LIA acquire a greatly increased lifetime when BSE 
in used. In the simulation by Samuels 12j, the non-local 
interaction of vortex rings in laminar flow of He II results 
in overall parabolic superfluid velocity profile matching 
the normal-fluid velocity profile. In the problem of su- 
perfluid turbulence decay at length scales greater than 
the interline spacing, the non-local interactions of vortex 
lines result in Kolmogorov cascade of eddies mimicking 
the behavior of the classical turbulence even without the 
normal component [TsL IT^ [l5| . 

In theoretical studies of turbulent vortex dynamics, an 
analytic solution is normally not available, and numeric 
simulations are of crucial importance. Clearly, the com- 
putational price of a direct simulation of BSE is macro- 
scopically larger than that of LIA, since to perform one 
step in time for a given element of a line one has to ad- 
dress all the line elements. In this Letter, we observe that 
this extensive addressing is actually extremely excessive 
and can be easily eliminated at the expense of a control- 
lable systematic error. The physical reason for that is as 
follows. The interactions at a distance r are significant 
only for the Kelvin waves of the wavelength A ^ r. The 
characteristic time of evolution of a harmonic of the wave- 
length A scales like A^ , which means that during the time 
interval much shorter than A^ the configuration of these 
harmonics remains essentially unchanged. Correspond- 
ingly, the maximal reasonable frequency of addressing 
the distances ~ r from a given element of the vortex line 
scales like l/r^. This leads to a numeric scheme where 
the long-range interactions of Kelvin waves do not reduce 
the simulation efficiency. We employ this scheme to sim- 
ulate the Kelvin wave cascade on a single vortex line at 
zero temperature. Our data for the cascade spectrum are 
in an excellent agreement with the predictions of analytic 
theory 

Scale separation scheme. — To be specific, we consider 
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the zero-temperature dynamics of Kelvin waves on a 
straight vortex hne. We separate different scales of dis- 
tance in the integral in Eq. Q by writing BSE in the 
form 

S^Ql+Q2 + Q3 + ...+Qn + ■■■ , (6) 

where 

n ^ [ (sq - s) X dsp 
Qn\s) = — / — r -7^ — • (7) 

Strictly speaking, Eqs. ©"O are meaningful only for 
n ^ 1, since at distances ^ oq from the vortex core 
the velocity field is model-dependent. With the same 
accuracy, one may replace in Eq. © the n ^ 1 non-local 
terms with a local term. (As we do it in our simulation 
scheme. — See below.) 
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FIG. 1; Comparing functions Q„{t) and Q„+i(t) for a fixed 
vortex line element. The data are obtained by simulating 
Kelvin wave cascade with full BSE. 

The idea behind the separation of scales of distance in 
Eq. jnjl is clearly seen in Fig. 1, where we compare two 
functions, Qnif) and Qn+i{t), in the regime of Kelvin 
wave cascade. The characteristic period of variation of 
the function Qn+i{t) is approximately 4 times larger than 
that of Qn- This is precisely what one would expect on 
the basis of qualitative analysis, since the typical wave- 
length of harmonics contributing to Qn is A„ = 2"ao, 
which leads to the estimate 

T„ oc (8) 

for the typical time of variation of Qn- The contributions 
to Qn from the harmonics with A <g; A„ are suppressed by 
the denominator in Eq. (jT)) in combination with a random 
interference in the numerator. 

The behavior of Q„'s suggests the following numeric 
procedure. Each Qn{t) is approximated by a piecewise 



polynomial function Qn{t) based on an appropriate mesh 
{4""*}' i ^ 1,2,3,..., and corresponding approximate 
values of ~ t^"'') obtained by integrating BSE with 
Qn{s,t) Qn{s,t). In the present study, we take the 
values of Q„ calculated at the time moments and 

(n) 

tj , and use the linear extrapolation to get the approxi- 
mation Qn{t) for Quit) at any t in the interval 

The key point then is to properly choose the spacing 
of the mesh in accordance with the estimate ©: 

4^,-4")cx7A^ (9) 

where 7 is a global (n-independent) small parameter con- 
trolling the accuracy of the scheme. (The scheme is 
asymptotically exact in the limit of 7 0.) With such a 

scaling of — tj"^ the long-range interactions do not 
form a bottleneck for the efficiency of the numeric pro- 
cess. Indeed, to numerically evaluate the integral ^ one 
has to perform oc A„ operations, but the time interval 
between two integrations scales like oc A^, which means 
that on average this integration costs only oc 1/A„ opera- 
tions per unit time. To control the accuracy, we perform 
a test simulation (at a smaller system size) and compare 
the results to the simulation of full BSE. At a reasonably 
small 7, we see no difference (up to a statistical noise, 
like the one noticeable in Fig. I^J in the evolution of the 
distribution of harmonics. 

Kelvin wave cascade. — The decay of superfluid turbu- 
lence at T — > leads to a cascade of Kelvin waves gener- 
ated by vortex line reconnections As it was argued 
by one of us 7], in the limit oi R/ao 00, when the dy- 
namics at length scales i? is well captured by LI A, in 
the low-energy part of the inertial range there takes place 
a specific regime when the energy flux in the wavenumber 
space is assisted by local self-crossings of fractionalized 
vortex lines. Corresponding spectrum of Kelvin waves 
(Kelvons) has the form 

nkOik-^/{\iiky , (10) 

where rik is the Kelvon occupation number (see below) 
and v is some dimensionless constant of order unity. The 
special role of local self-crossings follows from the inte- 
grability of Eq. ^ (the friction is absent at T = 0). — 
The existence of extra constants of motion excludes the 
possibility of the purely non-linear cascade. Vortex line 
self-crossings lift these constraints and push the cascade 
down to smaller wavelengths. 

Vinen conjectured jS] that — apart from the logarith- 
mic denominator — the spectrum (|10() can be derived just 
on the basis of dimensional considerations and is thus 
insensitive to the particular cascade scenario. If true, 
this statement would be of crucial importance, since at 
large enough wavenumbers the cascade should radically 
change its microscopic character. From dimensional esti- 
mates it is readily seen that with increasing wavenumber 
non-local interactions eventually become dominant, and 
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the cascade is driven by purely non-linear dynamics of 
Kelvons. According to Vinen, the change of the regime 
should not affect the Kelvon spectrum. Recent numeric 
simulations P, 0| seem to support this idea. However, 
an inconsistency can be seen in a dimensional analysis. 
Indeed, the spectrum Uk oc fc^"^ implies that the typical 
amplitude b\ of the Kelvon turbulence with the wave- 
length A is of the order of A. This leads to a strong non- 
linearity of BSE, which allows one to estimate the kinetic 
time Tkin as the dynamic time corresponding to bx ~ A. 
It follows then that Tkin(A) oc A^. With this estimate for 
the kinetic time, the spectrum Uk oc is inconsistent 
with the A-independent energy flux (see, e.g., Ref. ^ for 
the relation between spectrum, kinetic time, and energy 
flux). In our recent paper 0, we argue that the regime 
of the pure Kelvin wave cascade is characterized by the 
inequality 



aA = y < 1 , 



(11) 



with the parameter ax becoming progressively smaller 
with decreasing A. This small parameter allows one to 
develop a quantitative kinetic theory that predicts the 
spectrum 

Uk oc . (12) 

In terms of 6a, Eq. l(T^ is equivalent to 

6a oc A^/^ , (13) 

while the Vinen's spectrum means 6a oc A, that is aA 1 
for any A. 

Since the number 17/5 = 3.4 is rather close to 3, the 
accuracy of the data of Refs. P, turns out to be in- 
sufficient to distinguish between the two exponents. 

Numeric simulation. — The main goal of our simulation 
is to observe the exponent 17/5, and thus to verify the 
theory of Ref. ^. 

We consider a straight vortex line along the Cartesian 
z-direction with small (aA <^ 1) distortions in the {x,y)- 
plane; with the periodic boundary conditions. We use the 
z-coordinate as a parameter of the curve, and introduce 
the complex- valued function w{z, t) — x{z, t) -\- iy{z, t) to 
specify the configuration of the vortex line at a given time 
moment t. To proceed numerically, we employ the Hamil- 
tonian approach in terms of the discrete complex canon- 
ical variable Wj{t) = Xj{t) + iyj{t), j = 1,2,3, N, 
with the Hamiltonian function [mesh spacing is set equal 
to unity] 

AT 



N 



l + [(wj-+i - Wj-i)(w;%i - <_i)/8 + c.c] ^^^^ 



jj^i VU - 0^ + \wj -wi\ 

and corresponding equation of motion 

dH 



Here C is a constant of order unity. In the long-wave 
limit, the discreteness is irrelevant and the model 1)14(1 - 
(I15|l is equivalent to the Hamiltonian form of BSE in 
terms of the fmiction w{z,t), with z = j 0- From the 
theoretical point of view, the particular value of C, and 
the very presence of the first term in ((14|l . is not impor- 
tant. In practice, we introduce this local term to ren- 
der the model well-behaved at large momenta, where the 
spectrum of Kelvons is very sensitive to the short-range 
structure of the Hamiltonian. In our simulations, we used 
C = 23. 

The Kelvon occupation number is given by the square 
of the absolute value of the Fourier transform of Wj '5| : 



rik 



\Wk\ 



N 

Wk = N-'/'Y. 



Wi e 



ikj 



(16) 



(17) 



dw* 



(15) 



where k = 2'Km/N, m e [~{N - l)/2 , {N ~ l)/2] is an 
integer (we used odd N). 

There are different ways to generate a Kolmogorov- 
type cascade. Perhaps the most natural one is to cre- 
ate a steady-state regime with a source and drain, as 
it was done in the simulation of Ref. [J. However, an 
external source (random low-frequency force) is not nec- 
essary, since the cascade, by its definition, is a quasi- 
steady-state phenomenon, where the behavior of short- 
wavelength harmonics is adjusted to the energy flux com- 
ing from the long-wave harmonics. In a quasi-steady- 
state regime, the role of a source can be played by just 
the long-wavelength harmonics, as, e.g., in the simula- 
tion of Ref. llQ| . In both cases — with or without an ex- 
ternal force — the cascade regime requires some time to 
develop. We employ the following trick to circumvent the 
problem of transient period. We start with occupation 
numbers that progressively exceed the quasi-steady-state 
ones in the direction of shorter wavelengths. With such 
an initial condition, the cascade develops in the form of 
a wave in the momentum space propagating from the re- 
gion of higher wavenumbers towards the region of lower 
wavenumbers, the distribution behind the wave corre- 
sponding to the quasi-steady-state cascade 0|. For our 
purpose of resolving the spectrum (|12|l from the spectrum 
Uk oc k~^, it is convenient to use the latter distribution 
of occupation numbers as an initial condition. 

Introducing a drain — a decay mechanism for high- 
frequency harmonics — significantly improves the qual- 
ity of simulation, because the large-amplitude high- 
frequency modes distort the short-wavelength part of the 
cascade spectrum. We arrange the drain by periodically 
eliminating all the harmonics with \k\ > /ccutoff, with 
fccutoff = 7r/5. Note that the momentum fccutoff is 5 times 
smaller than the largest momentum in our system. We 
find it important to severely cut off the high-frequency 
harmonics because their dispersion strongly differs from 
the parabolic dispersion of Kelvin waves in the continu- 
ous space. 
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FIG. 2: A snapshot of time-averaged distribution of Kelvon 
occupation numbers. For time averaging we use time interval 
which is considerably shorter (by a factor of ~ 3) than the 
typical time of evolution. 

We simulated a system of the size N = 16385. Our 
main results are presented in Fig. |21 We do not show 
a large low-frequency part of the distribution with the 
spectrum Uk oc k~^. Having these ballast long- wave har- 
monics is important to guarantee that the inertial range 
is far away from the region where the wavelength is com- 



parable to the system size and the kinetics is strongly 
affected by the wavenumber discreteness. 

Fig. 121 reveals the wave converting the initial dis- 
tribution rife oc into a new power-law distribution 
Hk OC k~^. Quantitative analysis shows that the expo- 
nent (3 coincides with 17/5 with an absolute error less 
than 0.1, in an excellent agreement with the prediction 
of Ref. ^. 

Unfortunately, we were unable to check the prediction 
of Ref. for the relation between the energy flux and 
the amplitude A of the distribution nk = Ak~^'^/^. The 
problem is in calculating the mean value of the energy 
flux. With our quasi-steady-state setup, we can use only 
time averaging — with an essentially limited averaging pe- 
riod. This turns out to be insufficient to suppress huge 
fluctuations of the instant values of the energy flux. For 
these purposes, the steady-state setup of Ref. seems 
to be more adequate. 

In conclusion, we have proposed an efhcient numeric 
approach for simulating the dynamics of superfluid tur- 
bulence in the framework of Biot-Savart equation. The 
general principle is to use a less detailed addressing for 
longer-range interactions. The idea is implemented in a 
scheme for simulating Kelvin wave cascade on a single 
vortex line. With the new scheme we were able to ac- 
curately resolve the spectrum of the pure Kelvin wave 
cascade, observing an excellent agreement with the pre- 
diction of the kinetic theory of Ref. |2| . 
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